The distance duality relation from X-ray and SZ observations of clusters 
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duality relation between the angular and luminosity distances. We argue that such a test on large 
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I. INTRODUCTION 



Most cosmological observations provide compelling ev- 
idences that our universe is undergoing a late time accel- 
eration phase pj. However, there are still several debates 
about the physical interpretation of these observations. 
While it seems clear that the Friedmann equations for 
a universe only composed of normal matter (i.e. radia- 
tion and dust) even including dark matter cannot explain 
the current data, there are different ways of facing this 
fact. Either one can conclude that the interpretation of 
the cosmological data are not correct (i.e. we do not ac- 
cept the evidence for the acceleration of the universe, see 
Ref. for a recent critical review and e.g. Ref. 2]) or 
one tries to introduce new degrees of freedom in the cos- 
mological model. In this latter case, these extra degrees 
of freedom, often referred to as dark energy, can be in- 
troduced as a new kind of matter or as a new property 
of gravity. 

In the first approach one assumes that gravitation is 
described by general relativity while introducing new 
forms of gravitating components, beyond the standard 
model of particle physics, which must have some ef- 
fective negative pressure to explain the acceleration of 
the universe. Various candidates such as a cosmologi- 
cal constant, quintessence 0, Q with many potentials, 
K-essence |4| etc... have been proposed. But, one is 
still left with the cosmological constant problem |]| (why 
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is the density of vacuum energy expected from particle 
physics so small?) as well as the time coincidence prob- 
lem (why does the dark energy starts dominating today?) 
unsolved. From a cosmological point of view, these mod- 
els are characterized by their equation of state which can 
be reconstructed from the function E(a) — H 2 {a)/H^ 
where Hq is the Hubble constant at present and a the 
scale factor, either using the observation of background 
quantities or the growth of cosmic structures @ . 

The other route is to allow for modification of gravity. 
This means that the only long range force that cannot 
be screened is assumed to be not described by general 
relativity. Once such a possibility is considered, many 
classes of models exist (see e.g. Ref. @). For instance, a 
light scalar field can couple to matter leading to models 
of extended quintessence Q and more generally to scalar- 
tensor types of theories. Such theories have some difficul- 
ties to explain the current cosmological observations 
without a quintessence-like potential for the scalar field 
or a cosmological constant. This scalar field may also be 
at the origin of some variation of the fundamental con- 
stants, depending on its couplings, and violation of the 
universality of free fall (see Ref. [jjj] for a review). Other 
possibilities include brancworld models in which the stan- 
dard model fields are localized on a 3-dimensional brane 
embedded in a higher dimensional spacetime. Among 
braneworld models, a subclass of models have the prop- 
erty of allowing for deviation from 4-dimensional Ein- 
stein gravity on large scales. This is for example the 
case of some multi-brane models 0, multigravity [ l2 | . 
brane induced gravity or simulated gravity Q. In 
such models, gravity is not mediated only by massless 
gravitons, one therefore expects to have deviations from 
Newton inverse square law on large scales. Testing the 
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FIG. 1: Summary of the different classes of models and of the specific tests that can help distinguish between them (besides 
the equation of state and the growth of cosmic structures). The classes differ according to the kind of new fields and to the 
way they couple to the metric and to the standard matter fields. Upper-left class consists of models in which a new kind 
of gravitating matter is introduced, e.g. quintessence. In the upper-right class, a light field induces a long-range force so that 
gravity is not described by a spin-2 graviton only. This is the case of scalar-tensor theories of gravity. In this class, Einstein 
equations are modified and there may be a variation of the fundamental constants. The lower-right class corresponds to models 
in which there may exist massive gravitons, such as in some class of braneworld scenarios. These models predict a modification 
of the Poisson equation on large scales. In the last class (lower-left), the distance duality relation may be violated. 



Poisson equation on large scales may be a way to distin- 
guish between these alternatives |lall6|. 

The different types of models are summarized schemat- 
ically on Fig. 2] A diagnostic of the cause of the accel- 
eration of the universe will require to make many tests. 
In particular, the reconstruction of the function E(a) (or 
equivalently measuring the effective equation of state of 
the dark energy) will not be sufficient to distinguish be- 
tween many models. It is thus important to simulta- 
neously check for the Poisson equation, the growth of 
structure and the variation of the constants. 

Among these tests it has recently been pointed out in 
Ref. that the reciprocity relation and the distance du- 
ality relation that derives from it have also to be checked. 
The reciprocity relation is a relation between the source 
angular distance, r s , and the observer area distance, r a . 
The former is defined by considering a bundle of null 
geodesies diverging from the source and which subtends 
a solid angle dO s (see Fig- EJ). This bundle has a cross 
section d5 s and the source angular distance is defined by 
the relation 



The observer area distance r Q is defined analogously by 
considering a null geodesic bundle converging at the ob- 
server by 

dS* Q = r 2 Q dn o . (2) 

It can be shown that if photons travel along null geodesies 
and the geodesic deviation equation holds then these 
two distances are related by the reciprocity relation (see 
Ref. for a derivation) 

rl=rl(l + zY, (3) 

regardless of the metric and matter content of the space- 
time. Unfortunately, the solid angle df2 s cannot be mea- 
sured so that r s is not an observable quantity. But, it can 
be shown that, if the number of photons is conserved, the 
source angular distance is related to the luminosity dis- 
tance, Dl, by the relation 

D L = r s {l + z). (4) 

It follows that there exist a distance duality relation 



dS s = r 2 dQ s . 



(1) 



D L =D A {l + z) 2 



(5) 
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FIG. 2: A bundle of null geodesies diverging from the source 
(dash) subtending a solid angle dS7 s has a cross section dS s at 
the observer while a bundle converging at the observer (plain) 
subtending a solid angle df2 has a cross section dS at the 
source. 



that holds between the angular distance Da, the lumi- 
nosity distance Dl and the rcdshift z. This relation can 
be checked observationally . 

While the reciprocity relation holds as soon as photons 
follow null geodesic and that the geodesic deviation equa- 
tion is valid, the distance duality relation will hold if the 
reciprocity relation is valid and the number of photon is 
conserved. In fact, one can show that in a metric theory 
of gravitation, if Maxwell equations are valid, then both 
the reciprocity relation and the area law are satisfied and 
so is the distance duality relation (see Ref. [l8|l. 

There are many possibilities for one of these conditions 
to be violated. For instance the non-conservation of the 
number of photons can arise from absorption by dust, but 
more exotic models involving photon-axion oscillation in 
an external magnetic field [l!j can also be a source of 
violation [2(|. Note also that in principle both the reci- 
procity and distance duality relations hold for infinitesi- 
mal light bundles so that gravitational lensing may be a 
source of violation for macroscopic bodies. More drastic 
violations would arise from theories in which gravity is 
not described by a metric theory and in which photons 
do not follow null geodesic. 

In this paper, we propose and explore a potential new 
test of the distance duality relation based on Sunyaev- 
Zel'dovich |2l| (SZ) and X-ray measurements of clusters 
of galaxies. In Section [nj we first show that when the 
relation JSJ does not hold, cluster data do not give a 
measurement of the angular distance but of Da/ r] 2 where 



(6) 



In Section ITTTl we use existing cluster data to search for 
any hint that r\ = 1 may be excluded and then we discuss 
in Section IIVI the possibility to improve the accuracy of 
the test. 



II. RECIPROCITY RELATION FROM GALAXY 
CLUSTER OBSERVATIONS 

Galaxy clusters are known as the largest gravitation- 
ally bound systems in the universe. They contain large 
quantities of hot and ionized gas which temperatures are 
typically 10 7-8 K. The spectral properties of intra-cluster 
gas show that it radiates through bremsstrahlung in the 
X-ray domain. Therefore, this gas can modify the Cos- 
mic Microwave Background (CMB) spectral energy dis- 
tribution through inverse Compton interaction of pho- 
tons with free electrons. This is the so-called SZ effect. 
It induces a decrement in the CMB brightness at low 
frequencies and an increment at high frequencies. 

The possibility of using the SZ effect together with 
X-ray emission of galaxy clusters to measure angular dis- 
tances was suggested soon after the SZ effect was pointed 
out (see for example Ref. 22]). Used jointly, they provide 
an independent method to determine distance scales and 
thus to measure the value of the Hubble constant (e.g. 
Ref. mm for details). 

In brief, the method is based on the fact that the CMB 
temperature (i.e. brightness) decrement due to the SZ 
effect is given by 



AT S z ~ Ln e T e 



(7) 



where the bar refers to an average over the line of sight 
and L is the typical size of the line of sight in the cluster. 
T e is the electron temperature and n e the electron den- 
sity. Besides, the total X-ray surface brightness is given 
by 
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(8) 



where the volume V of the cluster is given in terms of its 
angular diameter by V 



D 2 A 9 2 L. It follows that 
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(9) 



The usual approach |22j , is to assume the distance dual- 
ity relation (rj = 1) so that forming the ratio AT| Z /Sx 
eliminates n e . Then, using a measurement of the angular 
diameter of the cluster and Eq. (|14|l one gets an estimate 
of the angular diameter distance and thus the Hubble 
constant. As a first conclusion, we point out that this 
method determines the angular distance only if the dis- 
tance duality relation is valid. Therefore one needs to be 
careful when using such data to test the distance duality 
relation. 

To make this point more precise, let us come back to 
the details of the method assuming the classical /3-modcl 
for the galaxy cluster [25| , that is assuming that the elec- 
tron density of the hot intra-cluster gas has a profile of 
the form 



n e (r) = n Q 



-i -3/3/2 



(10) 
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for < r < -Rciustcr and otherwise, i? c iuster being the 
maximum extension of the cluster. The temperature 
decrement due to the SZ effect in the Rayleigh- Jeans part 
of the spectrum is given by 



AT SZ (8) = -2-^V T 



(11) 



where we have assumed that the temperature of the hot 
gas, T e , is independent of r [T = T e (r = 0)]. 2£ max is 
the length of the path along the line of sight inside the 
halo of the cluster and 9 is the angular radial position 
projected on the celestial sphere from the cluster center. 
The X-ray emission is due to thermal bremsstrahlung and 
the surface brightness in a beam of angular diameter 59 
takes the form 



9 (9) D * 



dL 
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dV 



df 



(12) 



where the emissivity in the frequency band [1/1,1/2] is 
given by 



dL 



x 



dV 



(13) 



a(T e , V\, V2, z) is a function that depends on the proper- 
ties of the free-free emission for ions, on the mass frac- 
tion of hydrogen and on the gas temperature (see e.g. 
Refs [23, |24( f° r it s expression) . Introducing the angle 9 C 

by 



(14) 



where r c is the cluster core radius, and using the profile 
(|10|l we obtain in the limit i? c iustcr —> 00 
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and 
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where B is the Euler beta function. Using the definition 
of 77 from Eq. © , this latter expression rewrites as 



q <a\ S9 2 rf(z) 2 /6/3-1 1 



4tt (1 + z) 4 
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21 (l"6/3)/2 



(17) 



As expected, ATg Z /Sx eliminates n e and gives a mea- 
surement of the core radius r c from which we can deduce 
the angular diameter distance through Eq. (|14|l . When 
77 7^ 1, what is thus extracted from the data is an estimate 
of r c = r c /if. 

It follows from this analysis that, if we do not assume 
the distance duality relation to hold, what is in fact deter- 
mined is D^ ata (z) = r c /9 c which differs from the angular 
distance. We thus have access to 



Dr*{z) = D A {z)M\z) 



(18) 



which reduces to the angular diameter distance only 
when the distance duality relation holds. 



III. METHOD AND DATA ANALYSIS 

Our method is straightforward once we have made 
the previous remark. Using a data set of angular dis- 
tances determined from the combination of X-ray and SZ 
measurements, we have access to {z, Da(z)/t] 2 (z)}. To 
get 77 one needs to know the angular diameter distance. 
One possibility, and probably the most robust, is to es- 
timate it from its theoretical expression in a Friedmann- 
Lemaitre universe 



D A h (z) = I 



K 



dx 



'1/(1+*) x 2 E{x) 
where x — 1/(1 + z) and Jk is defined by 



f K {u) = 



(19) 



Ku sinh y'—Ku 
u > F=F^ 



K 



(20) 



respectively for K — (l — f2° nat — Q,° A ) jc? positive, 
null and negative. The function E 2 {x) = H(x)/Hq is 
explicitly given, for a A-CDM model by 



E 2 (x) = Cl« 



1 - n" 



(21) 



where fi^ at and are respectively the present matter 
and cosmological constant density parameters. 
We estimate rj(z) as 



r,{z) = JD^/Df^. 



(22) 



The error bars on this quantity will be estimated by a 
combination of the data error bars and of the la error 
bars on the cosmological parameters as obtained from 
Ref. IH 



< a t = 0.29±0.07 



0.73±0.05, 



= 0.73±0.04. 

(23) 

Note that in the case of a detection of 77 ^ 1 , the interpre- 
tation of the signal is not trivial since a varying equation 
of state may be undistinguishable from a violation of the 
distance duality relation. In such a case going back to 
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cluster 



redshift 



(Mpc) 



MS 1137.5 + 6625 0.784 

MS 0451.6 - 0305 0.550 

CI 0016 + 16 0.546 

RX J1347.5 - 1145 0.451 

Abell 370 0.374 

MS 1358.4 + 6245 0.327 

Abell 1995 0.322 

Abell 611 0.288 

Abell 697 0.282 

Abell 1835 0.252 

Abell 2261 0.224 

Abell 773 0.216 

Abell 2163 0.202 

Abell 520 0.202 

Abell 1689 0.183 

Abell 665 0.182 

Abell 2218 0.171 

Abell 1413 0.142 
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TABLE I: The 18 clusters of the Reese catalog with their red- 
shift used in our analysis. D^ ata refers to the angular distance 
determined in Ref. |27|l assuming that the distance duality re- 
lation holds. It leads, once this hypothesis is relaxed, to a 
measurement of rj with la error bars. 
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FIG. 3: rj as a function of the redshift for the 18 clusters 
of the Reese et al. [27]| catalog. The error bars include the 
observational error bars as determined by Reese et al. and 
the uncertainties in the cosmological parameters. 



SNIa data to get Dl and of X-ray data to get Da/t) 2 
may help to break this degeneracy. 

We use the catalog by Reese et al. [27,] that contains 
18 galaxy clusters, with redshifts ranging from 0.142 to 
0.784, all observed in X-ray and SZ (see Table UJ. Com- 
bining these data as discussed in Sect. ITU together with 
the theoretical estimate of the angular distance, we get 
a measurement of rj(z) for each cluster, using Eq. (|22|) . 
The result is summarized on Table [TJ and is depicted on 
Figure 



The question is then whether this data set is com- 
patible with rj = 1 or not. As can be seen from the 
original data (see Ref. |27j ), the error bars on D^ ata 
are not symmetric. To derive the distribution of rj 
we proceed in the following way. We first assume 
that the data points are independent so that the like- 
lyhood L — p(^data _ _ ^dataj^-j can f ac t or i Z ed as 

L — Y{ i Pi{vif ata \ri). We then introduce S defined as 



S 



-21nL. 



(24) 



If the probabilities Pi are Gaussian, S reduces to the stan- 
dard x 2 - I n one dimension, a variation AS = 1 around 
the minimum of S will give the la error bar. To proceed, 
we need to know the probabilities Pi. Without any fur- 
ther information, we assume that they follow a Gaussian 

distribution that is that ?7i + ^ ! - corresponds to a proba- 
bility distribution function of the form 



P(x) 



p -(x-rn) /25 + 

t (25) 

„-(x- Vi ) /2S_ ^ 



X <Jji 



(see e.g. Ref. |2g). 

We perform this analysis using two data sets. The first 
set, labelled 1, contains all the clusters, while in the sec- 
ond, labelled 2, we have removed the point at z — 0.374 
that lies outside of the other data points. This point 
corresponds to Abell 370 that clearly shows an appar- 
ent bimodal shape in optical and X-ray data, making its 
modelling as a single spherical potential a likely over- 
simplification for our purpose. The result is summarized 
on the plot 0] where we have displayed the function S, its 
minimum and the la confidence level. Fig. [S] compares 
the probability distribution function of rj to a Gaussian 
distribution fitted to the data. We obtain from our anal- 
ysis that 



rj = 0.87. 



-0.04 
-0.03 



for the first data set and 



„ _ o qi+0 04 



(26) 



(27) 



for the second. 

Additionally we have analyzed, for the second data 
set (i.e. without Abell 370), separately the low redshift 
(z < 0.3) data and the high redshift data. We find for the 
low redshift set rj = 0.89_o'og and for the high redshift set 
rj = 0.95^0 07- It is noteworthy that the largest departure 
from rj = 1 is at small redshifts. These results and the 
one obtained for the whole data set suggest that there 
is no significant violation of the distance duality relation 
from combined X-ray and SZ measurements. 



IV. PERSPECTIVES AND CONCLUSIONS 

Testing for the distance duality relation and/or the 
reciprocity relation can give some insight on the puzzling 
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FIG. 4: The function S as a function of rj for the two data 
sets (1 in plain and 2 in dash). The vertical bars indicates for 
each set the position of the minimum and the la confidence 
interval defined by AS = 1. 
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FIG. 5: Distribution of r\ obtained from the data of Fig. [3] 
compared with a Gaussian fit (dash) (left=set 1 and rights 
set 2). 



apparent acceleration of the universe derived from cos- 
mological observations. To distinguish between various 
models, one needs several complementary tests to the re- 
construction of the Hubble parameter as a function of the 
scale factor (or equivalently of the equation of state). Ex- 
amples of such tests are the test of the Poisson equation 
on large scales, the test of the constancy of the funda- 
mental constants and the test of the reciprocity relation. 
Figure ^ illustrates how the combination of these tests 
can help in identifying class of models. 

We have then shown that observations of galaxy clus- 
ters offer a test of the distance duality relation. In par- 
ticular, using SZ and X-ray measurements of the same 
clusters give an estimate of Da/tj 2 - An important con- 
sequence is that X-ray/SZ combined analysis does not 
give a measurement of the angular distances when the 
distance duality relation is violated. 

Testing the distance duality relation was already pro- 
posed in Ref. |l7J. In that work, different sets of data 
were used such as type la supernovae data to get the 
luminosity distance, and the FRIIb radio galaxies, X- 
ray clusters and compact radio data to derive angular 
distances. Interestingly, a general three parameter form 
of rj(z) was proposed in Ref. 0, based on general ar- 
guments about the violation of the conservation of the 
number of photons. Using their data, [l7| found a 2a vi- 
olation of the distance duality relation, mainly caused by 



an excess brightening of SNIa at redshift larger than 0.5. 
This analysis also allowed to put constraints on system- 
atic effects, such as SNIa extinction or evolution, that 
may bias apparent magnitudes. 

The analysis of the Reese et al. [27| cluster catalogue 
has shown that rj — 1 is marginally consistent with the 
data. In our study, we have not searched for a fit of a 
general expression for rj(z). Our main concern was first 
to use X-ray and SZ combined measurements of galaxy 
clusters to check whether the critical value rj — 1 was 
compatible with the data. Although we found that a 
value of rj sightly lower than 1 is favored, drawing any 
conclusion on the possible discrepancy between the dis- 
tances as predicted in the concordance model and those 
determined by our X-ray/SZ combined analysis is pre- 
mature. There are indeed several systematic effects that 
may bias our derivation of i](z), like over-simplification 
of cluster symmetry (substructures, tri-axiallity), or of 
their temperature and luminosity radial profiles. For ex- 
ample, the clusters that deviate most from the r](z) — 1 
line on Fig. 2 are those that clearly show bimodal struc- 
tures from X-ray, SZ and optical images (Abell 370, Abell 
773 and Abell 1689). In contrast, those showing a sin- 
gle emission region with spherical shape lie very close to 
i](z) — 1. It is interesting to note that when the three 
most bimodal clusters are removed from the sample, the 
Reese et al. [27| remaining clusters lead to 

»? = 0.93t°; ot (28) 

which is compatible to rj = 1 at a 2a level. The shape, 
temperature distribution etc., are key points to control 
in order to reduce the systematics that limit the accuracy 
of this test. 

Therefore, we will have to make sure that the marginal 
trend 1 77(2) < 1 survives further explorations of this 
method with new data. This trend is indeed related to 
the fact that X-ray/SZ analysis systematically favors a 
rather low value of the Hubble constant. Besides, the 
analysis of the low and high redshift subsets and the fact 
that the largest departure from 77 = 1 is at low redshift 
suggest that there is no violation of the distance dual- 
ity relation. In particular, parametric forms, such as the 
one proposed in Ref. 01 > predict a cumulative effect with 
redshift. 

More specifically, a larger number of clusters spread 
over the whole redshift range and showing simple appar- 
ent geometry (i.e. as compact and spherical as possible) 
must be selected carefully. It will improve to lower sys- 
tematics, to reduce errors bars on cluster data and, in 
turn, to provide much better angular distance estimates, 
making the test of the distance duality relation from X- 
ray and SZ measurement an efficient method. 



Note als o that this trend is opposite to the one obtained in 
Ref. |17|I (see their Fig. 1), which strengthens that there is in 
fact no systematic trend toward rj < 1. 
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